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The use of the reciprocal theorem has been shown to be a powerful tool to obtain the swimming 
velocity of bodies at low Reynolds number. The use of this method for lower-dimensional swimmers, 
such as cylinders and sheets, is more problematic because of the undefined or ill-posed resistance 
problems that arise in the rigid-body translation of these shapes. Here we show that this issue 
can be simply circumvented and give concise formulas obtained via the reciprocal theorem for the 
self-propelled motion of deforming two-dimensional bodies. We also discuss the connection between 
these formulae and Faxen’s laws. 


I. INTRODUCTION 

The locomotion of microorganisms in viscous fluids is constrained by the subdominance of inertial forces compared 
to viscous forces in the flows they generate. Some swimming strategies, such as reciprocal motion (in other words, those 
which display a time-reversal symmetry) can be effective at large scales while largely fruitless at small scales due to 
the kinematic reversibility of the Stokes equations [1]. In order to break time-reversal symmetry, microorganisms often 
pass deformation waves [2]. Determining the swimming kinematics of an organism undergoing a periodic deformation 
has been a subject of study since Taylor’s classic work on the swimming speed of an inhnite, two-dimensional sheet 
passing waves of transverse displacement [3], and continues to be an active area (see for example two recent reviews 
[4, 5], and references therein). 

Determining the motion of a swimmer in a fluid can be substantially simplified by appealing to the Lorentz reciprocal 
theorem [6]. Stone and Samuel showed that the swimming kinematics may be found without solution of the flow field 
generated by the swimmer if one has knowledge of the stress on a body of the same instantaneous shape undergoing 
rigid-body motion only [7]. This tool has been subsequently used to simplify a large number of problems from jet 
propulsion without inertia [8], to electrokinetic flow enhancement [9], to the motion of surfactant covered droplets in 
a background flow [10], as a few examples amongst many. 

Like Taylor’s swimming sheet calculation, many models use simplified geometries or lower dimensional geometries 
to gain physical understanding or mathematical tractability (for example [11-13]). The use of the reciprocal theorem 
for swimming, following Stone and Samuel [7], requires the solution of the stress under rigid-body motion which poses 
a problem for lower dimensional geometries where rigid-body motion is either ill-posed or undefined. However, unlike 
the case of rigid-body motion, force-free swimming is well dehned even for lower dimensional geometries, and we show 
here that these problems may be circumvented within the reciprocal theorem framework as well. We also show that by 
applying the reciprocal theorem to surfaces which are not material we can easily obtain the swimming kinematics for 
bodies whose shape is given by small deformations from a simple reference surface. Finally we rederive the formulas 
obtained for swimming via a direct integration of boundary integral formulation of the Stokes equation and discuss 
the connection with Faxen’s laws. 


II. SWIMMER MOTION 
A. Boundary motion 

In order to swim, a body undergoes a periodic deformation of its surface S{t). To describe the swimming gait of 
the body we use a virtual [14] (or unlocated [15]) body which does not move in response to the fluid. A point on the 
surface of the virtual swimmer, S{t), is given by fs(t) and the deformation of this surface in time is referred to as the 
swimming gait. Periodic deformations may be described as deviations from a reference surface 

^^(ro,^) - fo = Z\f(fo,t), (1) 

where fo is a point on the reference surface, Sq, of the unlocated body. When such a deforming body is submerged in 
a fluid, stresses exerted by the fluid lead generally to a rigid-body translation and rotation of the body. The position 
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of a point on the surface of the (located) real swimmer, X 5 S S{t)^ is hence given by the translation, by Xc(t) and 
rotation by ©(t) of that position on the virtual body, fg (see Fig. 1), hence 

xs(fo,t) = Xc(t) + 0 • fs(fo,t)- ( 2 ) 

A point on the reference surface in the lab frame, xq € So, is similarly writen Xq = Xc + © • fo- For simplicity we 
introduce the notation rg = © • fg and Tq = © • fo- 

Upon differentiation of the position of a point on the body we obtain the velocity 

u(xs) = =lJ + ftxrs + u^. (3) 

We have dehned U = and u'® = © • while the rotation operator is related to the angular velocity, fJ, by 
the relationship =0x0. The first two terms represent rigid-body motion while the third term represents 
the deformation of the surface due to the swimming gait [16, 17]. We note that the instantaneous translation and 
rotation of the body depends on the choice of reference point on the body, a satisfying resolution of this ambiguity is 
to specify that Xc be the center-of-mass of the body and then require that the gait of the virtual swimmer conserve 
momentum [14], but this complication is often unnecessary, in particular if one is ultimately interested in simply the 
time-averaged swimming speed. An alternative and equivalent description is to use a body-fixed basis to describe the 
deformation of the swimmer [18]. 

We describe here the locomotion of microorganisms small enough such that the Reynolds number of the flows 
generated may be taken to be zero and hence the Stokes equations are applicable 

V • cr = — Vp - 1 - TyV^u = 0, V • u = 0, (4) 

and the bodies are instantaneously force and torque free, 

F = / n-crdS' = 0, (5) 

Js 

L= f rs X {n ■ a) dS = 0, (6) 

JS 



FIG. 1. Schematic representation of a general swimmer. A position on the surface, S, of a swimmer in the lab frame is xs 
which is found by a translation (by Xc) and rotation of a point rs on a virtual swimmer. The instantaneous velocity of a point 
on the surface is given by translation U, rotation ft x rs, and deformation u®. 






3 


where the surface S' is a function of time and the normal to the surface n points into the fluid while cr is the stress 
tensor. For compactness we use six dimensional vectors which contain both force and torque F = [F L]^ and 
translation and rotation U = [U 

Due to the linearity of the Stokes equations for a Newtonian fluid, the flow field arising from a force- and torque-free 
body undergoing a surface deformation may be decomposed into two parts u = -|- u^. One part in which the body 

is undergoing only rigid-body motion U£i(x 5 ) = U -|- O x rg, and another part in which the boundary motion is due 
to the swimming gait alone, U 7 ’(x 5 ) = u"® [17]. There is, generally, a net force and torque in each problem. Under 
rigid-body motion, the drag Fo = —R U, where R is the rigid-body resistance tensor, while if the deforming body is 
held fixed {U = o), it generates an instantaneous ‘thrust’ (or swim force [19]), Ft- The sum of these forces must be 
zero, Fjy + Ft = o and so we write 

(;=/?-!. Ft. (7) 

In order to utilize this formula one must solve both the rigid-body motion problem, Uti (to obtain /?), and the thrust 
problem, ut (to obtain Ft). However, in a landmark paper [7], Stone and Samuel showed that by appealing to the 
reciprocal theorem [6], the solution of the flow field ut may be avoided entirely. 


III. THE RECIPROCAL THEOREM 

We denote ctt as the stress tensor associated with the velocity field due to rigid-body motion, Ut, while ctt as the 
stress tensor for associated with the velocity field due to deformation, ut- The stress and velocity fields satisfy the 
relationship V • ctt • ut = V • ctt • ut = 0 and by integrating over the fluid volume exterior to S by means of the 
divergence theorem [7], one finds 


Ft U= [ n-(TD-u^dS. (8) 

Js 

Due to the linearity of the Stokes equations we may write cro = T ■ U and hence the resistance tensor is given by the 
relationship 


R = - 


n T 

rs X (n • T) 


dS*. 


Substituting into (8) we find 


Ft = [ • (n • T) dS', 

Js 


(9) 


( 10 ) 


indicating that the hydrodynamic ‘thrust’ Ft may be written as an integral of the rigid-body stress tensor weighted 
by the gait u"®. Finally substitution into (7) gives the rigid-body motion of the swimmer 


U = R~^ 



{n-T)dS , 


( 11 ) 


as shown by Stone and Samuel [7]. In this form we see that the swimming velocity is a simple linear functional of 
the gait, U = >C(u'^), from which well known constraints on swimming in Stokes flows such as the scallop theorem or 
the rate invariance of the distance travelled in a stroke follow directly [14]. Note that for two-dimensional bodies in 
Stokes flows the resistance, or dragging, problem may be ill-posed; the stress on the surface of a translating cylinder, 
for example, is singular in the Reynolds number, nevertheless the product n • T • R^^ is well defined and hence this 
formalism may be used for lower dimensional bodies. Finally we may recast the swimming motion with respect to 
the unoriented body 


u = eu = eiR~^- 


drs 

dt 


(n • f ) dS' 


where U = [U O]^ indicates unoriented rigid-body motion and the six dimensional rotation operator. 


0 = 


© 

0 


0 

© 


( 12 ) 


(13) 
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IV. SIMPLE BODIES 

Analytical progress can be made by considering shapes which align with a coordinate system. We describe below 
the swimming of bodies with spherical, cylindrical and planar shape. The spherical case was demonstrated by Stone 
and Samuel [7]. The cylindrical and spherical cases show that the undefined or singular difficulties with the rigid-body 
resistance problem can be circumvented to yield formulas equivalent to the spherical case. 


A. Spherical swimmers 

If the body is a sphere of radius a, then the rigid-body problem is well known; the resistance matrix is diagonal 
and hence easily invertible while n • T = — |y[I 2H] where H = —r x I where I is the identity tensor. With these 
formulas the swimming velocity for a sphere with a prescribed velocity field, u'®, on its surface is 

^ 47ra2 Jg 

(where the angle brackets denote an average over the surface S) as shown by Stone and Samuel [7]. 


I 

3 ctT 

2a^ “ 


U'^ dS” : 


3 

2a2 


(rs X u'^) 


(14) 


B. Cylindrical swimmers 


Consider the flow around a rigidly translating and rotating cylinder of radius a. The cylinder can translate in any 
direction but we restrict rotations to the axial direction (r2j[). The surface traction may be written as 


„.T = -^[/3I|l+aIu 2H.I„] 


(15) 


We use the || subscript to denote the axial direction while T denotes components perpendicular. For translation 
normal to the surface of the cylinder, a is a dimensionless constant which is singular in Re [20]. For axial translation 
we may say we have an outer cylinder of radius R at infinity and so j3 = limfl_>oo Integrating over the surface 

of the cylinder we obtain the resistances per unit length 


p-1 ^ ^ T 

47r77a2 

Combining these terms we obtain the swimming velocity 


U = - 



(nS) 

(rs X u^) ■ I|| 


(16) 

(17) 


(18) 


Notice that the result is independent of a and /3 and thus the singular nature of the resistance problem is avoided 
because the product n • T • is well defined. Squires and Bazant first derived the form of the solution for the 
translation of an infinite cylinder perpendicular to its axis, in the context of electrophoretic migration, by noting that 
n ■ (Tg) is constant [21]. 


C. Planar swimmers 

Consider a flat two dimensional sheet which may have a different prescribed velocity on each side of the sheet, 
unequally spaced between two rigid walls. If the sheet can only move in its own plane Uy, the resistance problem 
would be shear flow where 

n.<Tp = -^U||. (19) 

If the distances between the two walls are hi and ft .2 we have then as the only non-zero mobility 


R-fu “ [^1 ^ + ^2 I|h 


( 20 ) 
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and with this Eq. (11) becomes 


U = 


-1 


hi + /i2 

In an unbounded fluid, /ii = ^2 —>■ oo, we obtain 

1 ^ 


^2 (u‘^^) +/ll (u'^2) y 


U = - 


.((u^0 + (u"^>)=- ^1' •(«">. 


( 21 ) 


( 22 ) 


If the the sheet is symmetric, and one need only solve the half space problem. 


D. General form 


All the formulas derived above for finding the swimming velocity of an unbounded sphere, cylinder or sheet, U due 
to a prescribed boundary condition, u'^, take the general form 


C (rs X u 


">J 


U = - 

where C is simply a geometric factor. This may equivalently be written as 


(23) 


(u) = 0, 
(rs X u) = 0, 


(24) 

(25) 


where we see that because of the symmetry of the stress, both the surface-averaged velocity and the surface-averaged 
first moment of the velocity must be zero. This connection was demonstrated for a sphere by Felderhof [16]. Equations 
(24) and (25) may be obtained through a direct integration of the boundary integral formulation of the Stokes equations 
over the surface of a sphere, or a cylinder following Batchelor’s approach to the derivation of Faxen first law (see 
Appendix) [22]. For a flat sheet, we observe that the force on the surface is given simply by 


F = r] 



= 0 , 


(26) 


which is indeed true on every plane parallel to the sheet and so we must have 


(u) = 0, 


(27) 


on the surface of the sheet. The constant of integration is determined by the far-field condition. 

Recall that u(xs) = U -|- x rg -|- u'® and so upon substitution into (24) and (25) we find 

(28) 
(29) 


U = -(u^), 

(I — nn) ■ ft = — (rs x u'^) . 


For a sphere (I — nn) = |l and so C = 3/2a^ in Eq. (23), while for a cylinder where H = riy, hence (I — nn) FI = fl 
and so C = 1/a^. 


E. Deforming bodies 

We have described above how to obtain the motion for objects with simple shapes which align with coordinate axes. 
If the time dependent surface of a swimmer, S(t), deviates only slightly from a simple reference surface, Sq, then we 
can, through Taylor series recast the problem to obtain the velocity field on Sq [23, 24]. The velocity field is expanded 

u(xs) = u(xo) -f (xs - Xo) • Vu]^^ -f 0{\At\^), (30) 

then rearranging we obtain the boundary condition on S'o, 

u(xo) = U -I- H X ro -I- u‘®“. 


(31) 
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where 

f) /\\' 

+ n X Z\r - Z\r • VuL + 0(\Ar\^), (32) 

ot ° 

and Z\r = rg — rp. We remove the dependence on orientation of this instantaneous boundary condition by taking the 
product of both sides with to obtain 

- ^ y\ -p 

u(ro) = = ^—+ X Z\f - Z\f • Vu +C>(|Z\f|2), (33) 

at Vo 

where Cl — ©^ • fl and u(f) = ©^ • u(x) is the velocity field with respect to the unlocated body with x = + © • f . 

This formulation leads to a swimming problem defined entirely on Sq which satisfies the correct boundary conditions 
on S. In order to solve for the swimming kinematics, we apply the reciprocal theorem on Sq, namely we decompose 
the problem into finding the resistance to the rigid-body translation and rotation of Sq and the thrust force generated 
by u'^°. Nothing in the formulation of the reciprocal theorem for swimming prevents its application on a surface that 
is not material and by doing so one avoids the expansion of a surface integral over a non-trival domain [25] and in 
addition now we need only to solve for a single resistance problem. For bodies that are nearly spherical, cylindrical 
or planar in particular we may write 





1 

o 

_1 

u = - 

_C(rg X u^®)_ 

= -©• 

^(Jfg X 


= e0. 


(34) 


One may first determine U and then integrate fl to obtain the orientation in the lab frame. Note that the swimming 
gait on So, u'®®, depends on gradients of the (unknown) flow field u and the rotation rate of the swimmer 17. However 
if we take Z\f = X]m>o where e <C 1 is a dimensionless measure of gait amplitude and then consequently expand 
the velocity field, u = J2m we overcome these difficulties perturbatively, at each order 




u 


Uo 


dr I 

Ik’ 
dr 2 ^ 

— -h Oi X ri - ri • Vui 


(35) 

(36) 


This yields the instantaneous boundary condition for the velocity field on Sq in the orientation order by order and 
then by Eq. (34) we obtain U. 


V. EXAMPLE MODEL SWIMMERS 

In the following section we revisit two classic model swimmers, Taylor’s swimming sheet [3] and a squirming [26] 
cylinder [27], using the reciprocal theorem formalism presented above. 


A. Deforming sheet 

We describe here an infinite two-dimensional sheet with time varying deformation of its surface from a planar 
reference surface. Due to its infinite extent the body does not rotate, ©(t) = ©(0), and we take ©(0) = I for 
convenience. If the deformation is in the form of periodic waves then we may write 

Zir = A^c„e“^ (37) 

n 

The vector of complex Fourier coefficients, c„, delineates wave shape (cg = 0), while the wave variable, z = k-rg —wt, 
depends on the wave vector, k, and the actuation frequency, uj. The boundary conditions for an extensible sinusoidal 
sheet. Cl = —ijl&y, and an inextensible sheet were first described by Taylor [3], and catalogued in detail elsewhere 
[28]. 

Non-dimensionalizing lengths by jkj and time by uj we obtain 

Z\r'= e^c„e“^ = er'i, (38) 
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where e = ^|k| and primes indicate dimensionless quantities. Dropping the primes and assuming dimensionless 
variables we obtain the boundary conditions on Sq order by order in e, 

uf“ = 

n 

n 

We note that for this general deformation, the values of u'^‘> must be computed for each side of the sheet S'! and S '2 
as the value of the mean need not be equal. 

Expressing the swimming velocity U(e) = and using Eq. (22) for an unbounded swimmer we see imme¬ 

diately that due to the wave-like deformation 

Ui = =0. (41) 

To find the swimming speed at second order we must obtain the complete first order flow field in order to evaluate 
the gradients Vui. The exception is if the sheet is only deforming tangentially in the plane of Sq, in that case only 
in-plane derivatives, Vq, are required 


(39) 

(40) 


Voui = Voi-i = n2kc„e“T (42) 

n 

The reciprocal theorem for an unbounded swimmer, Eq. (22), then yields 

U 2 = - • k (43) 

n 

where ^ denotes the complex conjugate. Interesting to note is that when the sheet is subject to tangential deformations 
only, the presence of walls is strictly irrelevant, provided the prescribed velocity field is the same on both sides of the 
sheet, as described in Eq. (21). 

To include out-of-plane deformations of the sheet we must solve for Ui. At this point we’ll assume only fully two- 
dimensional flow fields with c„ = an^x -b bn^y where e^, = k/|k| while By is normal to the sheet. The flow field may 
be obtained by means of the streamfunction u = and the general case follows directly from classical examples 

for single modes as each Eourier mode is decoupled at leading order (details are hence omitted). The mean of the 
boundary condition on each side of the sheet is given by 

) = “ X! ^ A(n/i)ia„&j^ - B{nh)bnb'l] e^,, (44) 


where the sign corresponds to the surface with normal n = ±ey, and 


A{x) 


sinh(2a;) — 2x 
sinh^(a;) — ’ 


B{x) 


sinh^ (x) -b 
sinh^(x) — x^ 


(45) 


This difference in the mean value of the boundary condition on each side of an asymmetric sheet was recently noted 
by Felderhof [29]. Using Eq. (44) we obtain the swimming speed of a sheet that is unevenly spaced between two walls 
directly from Eq. (21) 


U2 = E^ 


a„aj, -b 


h2A{nhi) — hiA{nh2) . h2B{nhi) + hiB[nh2) 

hi + h2 " hi + h2 


bnbi 


G/ 7 .. 


(46) 


Without out-of-plane deformations, (6„ = 0), this result matches Eq. (43) and indeed the motion of the sheet is 
unaffected by the presence of a wall. It is important to note here that this result is inconsistent with that found by 
Katz [30] for a single mode swimmer near two walls due to the minus sign in the second term which arises from the 
asymmetric boundary condition noted in Eq. (44). Now with only out-of-plane deformations, a„ = 0, the sheet is 
symmetric and we obtain Reynolds result [31] for a single non-zero mode. In the particular case of an unbounded 
swimmer, h —>■ 00 , while A(nh) —)■ 2sgn(n) and B(nh) —)■ 1 and so the swimming velocity is given by 




n 


( 47 ) 
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This differs with the general result first found by Blake [32], again due to the asymmetry in the boundary condition 
on either side of the sheet, specifically because the term ±2isgn(n)a„6]j ultimately cancels out when we average the 
top and bottom of the sheet. Felderhof argues that for a finite circular disk the velocity difference across the disk 
induces a rotation of the body [29]. Taylor’s sinusoidal swimming sheet, where the only non-zero mode is bi = —il2, 
yields the classic result U = — to leading order [3]. 


B. Deforming cylinder 

Following Blake [27], we describe the deformation of the surface of a cylinder through material coordinates on the 
virtual swimmer rs{0o,t) = + eri{0o,t) and 9s{0Oit) = 0o + then a position on the material surface in 

this frame is simply is = rser{0s)- By expanding in e we obtain 

Z\f = ^ = e(rier -b ro6»iee) -b e^(ri6»iee - (l/2)ro0fer) + ... (48) 

m>0 


where the convention is = e^(^o)- Upon differentiation we have fi = riBr + T^Oieg and f 2 = {ri9i + riOi)eg — 
ro0i0ier. The motion of the body then follows from the equations (35) and (36), to leading order we have 




(rie^) -b To {0iee\ 

Tq ^ (yq X fi^ 


-1 

0) 


(49) 


A proper swimming gait is periodic in time and so a time-averaged velocity in the swimming frame is zero, Ui = o. 
Of course the time-averaged velocity of the swimmer need not be zero in the lab frame over the same period as 
rotation and translation are coupled. The leading-order rotation vanishes however, if the tangential deformation has 
zero mean, as we now assume. 

To obtain locomotion at quadratic order one needs to solve for the entire velocity field at leading order unless only 
tangential deformation occurs, in other words ri = 0, yielding a so-called squirmer. In this case we obtain 



(i2^ - ^fi • Vofi^ 


ro (0idea0iee'j 

r-2 

'o 

fo X fa - ^fo X (fi • Vofi)^ 




(50) 


Exploiting the periodicity of the geometry one can write the deformation 0i 


We see then that 


ro'^[bi{t)]ex + roRe[6i(t)]ej^ 

0 


while 


\ ' fi h 

2 Z^n 

+ (^i+l + 


inhnbie^ 


(51) 


(52) 


For directed motion it is typical to restrict 0i to be an odd function, doing so by taking = —{ij2)l3n where the 
coefficients /3„ = —/3_„ are real, one obtains up to quadratic order 


0 = 


e^e^, -b Y.n>0- Pn+l)^^: 

0 


(53) 


in agreement with Blake’s result [27]. Finally, the leading order time-averaged rigid-body motion of the swimmer in 
the lab frame is then given by 


U = 0(0) • 0 = 0(0) • 


e^TE„>o^/3n(/3n-i -/3„+i)e,^ 

0 


(54) 
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VI. CONCLUSION 

The reciprocal theorem for swimming allows the solution of the kinematics of a swimming microorganism without 
the resolution of the entire flow-held. Here we’ve shown the reciprocal theorem can also be applied to study the 
motion of two-dimensional swimmers. The reciprocal theorem shows that both the surface-averaged velocity and the 
surface-averaged first moment of the velocity must be zero on swimming sphere, cylinder and sheet, a result that may 
also be obtained directly by integrating the boundary integral formulation of the Stokes equations. 
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Appendix A: Boundary integral formulation (sphere) 

The reciprocal theorem was used to obtain the formulas 

(u) = 0, (rs X u) = 0. 

These formulas for a sphere and a cylinder may be alternatively found by using the boundary integral formulation of 
the Stokes equations following Batchelor’s derivation of Faxen’s first law [22]. We also include an arbitrary background 
flow for full generality. 


1. 3D boundary integral formulation 


With the boundary integral formulation of the Stokes equations for a point on the boundary X5 

u(xs) = Uoo(xs) - ^ [ [G(xs,y) • f(y)-f r 7 u(y) • T(xs,y) • n(y)] dS'(y) 

Airr] Js L J 

= U -I- O X rg -I- u'®, 

where Uqo is an arbitrary background flow, the Stokeslet, G, and its associated stress tensor, T, are given 

„ / , I XX „, , xxx 

where x = x — y. 

If we now average this velocity field over the surface S, written as (u), we obtain 

(u) = (uoo) - ^ (G) • f(y) + Wiy) ■ (T) • n(y) d5'(y) = V+ (u®) . 

Using two identities [22, 33], for averages over the surface of a sphere, 

1 


(T) . n(y) = ^I, 


we obtain 


U = (Uoo - u^) - 


1 


Girrja 


F. 


If the swimmer is force free, F = 0, moving in an otherwise quiescent fluid, Uqo = 0, then we have simply 

(u) = 0. 


(Al) 


(A2) 


(A3) 

(A4) 

(A5) 

(A6) 

(A7) 
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Similarly 



(rs X u) = (rs x Uoo) - f [ (n x G) • f(y) -b r 7 u(y) • (n x T) • n(y)] dS'(y), 

47r77 Js L J 

= ((rs • rs)I - rsrs) • + (rs x u'^) . 

(AS) 

Noting that rs = 

an and utilizing the following identities 



(n X G) = ^n(y) x I, 

(A9) 


(n X T) • n(y) = ^n(y) x I, 

(AlO) 

leads to 

>■"» 8„„3L = 1). 

(All) 

For a torque-free 

swimmer in a quiescent background flow we have simply 



(rs X u) = 0. 

(A12) 


The formulas above show that the rigid-body motion of a spherical swimmer is identical to that of an inert sphere 
placed into a background flow given by Uqo — u'^ according to Faxen’s laws. 


2. 2D boundary integral formulation (cylinder) 


With the boundary integral formulation of the Stokes equations for a point on the boundary xg, where S now refers 
to the one dimensional perimeter of the cylinder, we have 


u(xs) = Uoo - ^ [ [g(xs, y) • f(y) -b 77u(y) • T(xs, y) • n(y)j d5'(y) 


27r?7 


IS 


= U -b Jl X rg -b , 

where the 2D Stokeslet G and its associated stress tensor T are given 


T(x,y) = -4^ 


G(x,y) = -Iln(|x|) -b 


Averaging the velocity field over the surface of the cylinder using the identities 

(G)= [^-Ina] I, 


(T) • n(y) = -I, 


we obtain 


U = (Uoo - u"®) - 


s\ (5-lna) 


dTrry 


F. 


(A13) 

(AM) 

(A15) 

(A16) 

(A17) 


where F is a force per unit length. For a force free swimmer in an otherwise quiescent flow we have simply 

(u) = 0. (A18) 

Conversely taking the average {rs x u) with the identities 

(n X G) = n(y) x I, (A19) 


leads to 


(n X T) • n(y) = -n(y) x I, 


Si=4<rsx(u„-uq>-j^L. 


(A20) 

(A21) 


where L is torque per unit length. For a torque free swimmer in an otherwise quiescent field we have simply 

(rg X u) = 0. (A22) 
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